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I. INTRODUCTION 



Despite splendid successes of the standard model of particle physics, some fundamental 
questions remained unanswered, among which the origin of fermion masses is one of the 
most important problems. In the standard model, these masses are taken as free parameters, 
and can only be extracted from experimental data. Phenomenological analysis aiming at 
discovering the numerical relations between fermion masses is useful and practical for the 
exploration of the mystery of lepton and quark masses. There have been some conjectures 
on this problem (for example, [l|, |2|), in which Koide's relation is the most accurate 

one, which links the masses of charged leptons together, 
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m e + JmTa + Jm r ) 



where m e , m^, m T are the masses of electron, muon, and tau, respectively. In order to see 
the accuracy of Koide's relation, we can introduce a parameter ki, 

m e + m M + m T 



k,= 



With the data of PDG 



0.510998902 ± 0.000000021 MeV, 



0.000005 MeV and m T = 1776.99±g;l§ MeV, we can get the range of k h h 



(1) 



105.658357 ± 

_ -i +0.00002635 
~~ 1 -0.00002021? 



which is perfectly close to 1. 

This relation was speculated on the basis of a composite model of quarks and leptons 
and the extended technicolor-like model Q|. The fermion mass matrix in these models is 
taken as 

Mf = m f Q GO f G, 

where G = diag(gi,g 2 , #3). With the assumptions gi = g^+g\ ■ , Yl% di^ = an d ^2ii9i^) 2 — 
3(g^) 2 , and the charged lepton mass matrix is 3x3 unit matrix, Koide obtained this relation. 
This relation is deduced in several other different models by Koide also jSJ . (For a review , 
see 0j Several explanations for this excellent relation were presented in the past ten years. 
Foot [gj] gave a geometrical interpretation to it, 



cos 61 



1,1,1 



/m. 



\(y/m;,^/m^,y/m;)\\(l,l,l)\ V3^/m e + m^ + m T ' 
where Q\ is the angle between the points (y/m e , Jva^ y/m T ) and (1, 1, 1). Commonly, we can 



extend (1, 1, 1) to a more general mass square root vector (-v/mo, y/rriQ, v/^o), however, the 



value of mo does not affect explicitly the value of 61. Furthermore, we can find the relation 
between ki and 9i, 

k 1 



2 cos 2 61 



and 9i = f . 

Besides the miraculous success of Koide's relation for charged leptons, two further ques- 
tions emerge naturally, 

(1) . whether this relation can be applied to quarks and neutrinos, 

(2) . whether this relation holds well at different energy scales. 

Q 

The first question has been discussed by Esposito and Santorelli |9(, Rivero and 



Gsponer and us In Ref. we applied Koide's relation to quarks and neutrinos, 
and estimated the masses of neutrinos. For the second question, due to the renormalization 
effect, the masses of quarks and leptons vary with the energy scale of interaction, so we 
must examine Koide's relation at different energy scales explicitly, before one can take the 
extension of the relation to quarks and neutrinos seriously. 

In Section II, we check the application of Koide's relation to quarks and charged leptons. 
By using the previously introduced parameters k u , k^ to characterize the deviations of w-type 
and rf-type quarks from Koide's relation we find that k u and ka keep almost invariant 
in a very wide energy range. This means that Koide's relation is a universal result which 
is independent (or insensitive) of the running masses of quarks and leptons. In Section III, 
we apply Koide's relation to neutrinos, with both the normal and inverted mass schemes 
considered. Furthermore, with some analogies and conjectures between quarks and leptons, 
the neutrino masses are predicted. Finally, in Section IV, we give some discussions to Koide's 
relation. 



II. KOIDE'S RELATION FOR QUARKS AND CHARGED LEPTONS 



In order to see whether Koide's relation holds well at different energy scales, we must get 
the values of quark and lepton masses first. Because of the renormalization effect, the scale 



dependence of the running quark masses is determined by 



12, 
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where /i is an energy scale, n is the effective number of quark flavors 
effective QCD coupling constant, which is also a /i-dependent function, 
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, and a s is the 
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and Po, Pi, Pi are the coefficients of the renormalization group equa- 

fji— a s = P(a s ), 
djj, 
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Using Eq. (j2J), the numerical results of the running masses of quarks were obtained by 
Fusaoka and Koide [13], as summarized in Table 1, in which both low and high energy scales 
are calculated. 

First, for the energy scales lower than the spontaneous symmetry breaking energy scale 
Aw of the electroweak gauge symmetry SU(2)l ® U(1)y, seven different energy scales are 
taken into account, i.e., \i = 1 GeV, /1 = m c , /1 = m&, /x = my/, \i = mz, A* = m ti an d 
H = Ayr, where mw and mz are the mass of W and Z bosons. 

Second, for the energy scales extremely higher than Ay^, the evolution of the Yukawa 
coupling constants must be considered. In [l3j, two different models are presented. One 
is the standard model, and the other is the minimal SUSY model. In both of these two 
models, the Hamiltonian of the fermion fields can be written as 



H = Y^La^Raj^a + H.C., 



TABLE I: The masses of quarks at different energy scales. The upper seven rows in Table I are 
the cases of lower energy scales, the 8th and 9th rows are the cases of higher energy scales in the 
standard model, and the last two rows are the cases of higher energy scales in the minimal SUSY 
model. 



A* 


m u (MeV) 


m c (MeV) 


m t (GeV) 


m d (MeV) 


m s (MeV) 


m b (GeV) 


1 GeV 


4.88 


1506 


475 


9.81 


195.4 


7.18 


m c = 1.302 GeV 


4.18 


1302 


399 


8.40 


167.2 


6.12 


m b = 4.339 GeV 


3.17 


949 


272 


6.37 


126.8 


4.34 


m w = 80.33 GeV 


2.35 


684 


183 


4.73 


94.2 


3.03 


m z = 91.19 GeV 


2.33 


677 


181 


4.69 


93.4 


3.00 


m t = 170.8 GeV 


2.23 


646 


171 


4.49 


89.4 


2.85 


A w = 174.1 GeV 


2.23 


645 


171 


4.48 


89.3 


2.85 


10 9 GeV 


1.28 


371 


109 


2.60 


51.9 


1.51 


2 x 10 16 GeV 


0.94 


272 


84 


1.94 


38.7 


1.07 


10 9 GeV 


1.47 


427 


149 


2.28 


45.3 


1.60 


2 x 10 16 GeV 


1.04 


302 


129 


1.33 


26.5 


1.00 



where a — u,d, i,j = 1, 2, 3, Y°- are the Yukawa coupling constants, and 0° are the vacuum 
expectation values of the neutral components of the Higgs bosons (f> a . In the standard model, 
4>° = 4>u = 4>°di an d m the minimal SUSY model which has two Higgs bosons, 0° = 0° sin f3(fx) 
and 4>° d = 0° cos /?(//). The mass matrices at the energy scale fx = A w of quarks are given by 

m(ii) a = -^=Y({i) a v a , 

where v a are the vacuum expectation values of 0°, and v a = \/2(0°). In the standard model, 
v u — Vd = v^W, and in the minimal SUSY model, y/v \ + v \ = \^2A W . 

We can see from Table I that all the quark masses decrease with the increase of the energy 
scales, and the masses in the minimal SUSY model are a little bit larger than those in the 
standard model at the same energy scales. 

Now, we can examine whether Koide's relation holds well for quarks with the help of 
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TABLE II: 

x m Du an d k u at different energy scales. 





m u (MeV) 


m c (MeV) 


m t (GeV) 




Vu 




1 GeV 


4.88 


1506 


475 


313.75 


98958 


1.341 


1.302 GeV 


4.18 


1302 


399 


311.48 


95455 


1.338 


4.339 GeV 


3.17 


949 


272 


299.37 


85804 


1.333 


80.33 GeV 


2.35 


684 


183 


291.06 


77872 


1.328 


91.19 GeV 


2.33 


677 


181 


290.56 


77682 


1.328 


170.8 GeV 


2.23 


646 


171 


289.69 


76682 


1.327 


174.1 GeV 


2.23 


645 


171 


289.24 


76682 


1.327 


10 9 GeV 


1.28 


371 


109 


289.84 


85156 


1.335 


2 x 10 16 GeV 


0.94 


272 


84 


289.36 


89362 


1.339 


10 9 GeV 


1.47 


427 


149 


290.48 


101361 


1.347 


2 x 10 16 GeV 


1.04 


302 


129 


299.01 


124038 


1.359 



the data in Table I. To characterize the deviation of quark masses from the exact Koide's 



relation, here we introduce, as was done in Ref. 
in Eq. (P), 



two parameters k u and k d similarly as 



and 



k _ m u + m c + m t _ l+x u + y u 

f (y/rn^+ a/^+ y/m~t) 2 fll + v^+v^) 2 ' 

= m d + m s + m b _ l+x d + y d 

where x u = m c /m u , y u = m t /m u , x d = m s /m d , and y d = m b /m d . 

With the numerical results of quark masses in Table I, we can calculate all the x u , y u , 
x d and y d of both the u-type and d-type quarks at different energy scales, and then get k u 
and k d at different energy scales straightforwardly. These results are listed in Table II and 
Table III. 

We can see from Table II that both x u and y u decrease with the increase of fi at low 
energy scales and increase slightly at high energy scales. However, we find that in spite of 
the change of x u and y u , k u is almost invariant (approximately 1.33) at all energy scales, 
which means that k u is a constant independent of the running of quark masses. This is an 
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TABLE III: x^, yd and kd at different energy scales. 





m d (MeV) 


m s (MeV) 


m b (GeV) 


X<l 


Vd 


kd 


1 GeV 


9.81 


195.4 


7.18 


19.92 


731.91 


1.068 


1.302 GeV 


8.40 


167.2 


6.12 


19.90 


728.57 


1.067 


4.339 GeV 


6.37 


126.8 


4.34 


19.91 


681.32 


1.057 


80.33 GeV 


4.73 


94.2 


3.03 


19.92 


640.75 


1.048 


91.19 GeV 


4.69 


93.4 


3.00 


19.91 


639.66 


1.048 


170.8 GeV 


4.49 


89.4 


2.85 


19.91 


634.74 


1.046 


174.1 GeV 


4.48 


89.3 


2.85 


19.93 


636.16 


1.047 


10 9 GeV 


2.60 


51.9 


1.51 


19.96 


580.77 


1.032 


2 x 10 16 GeV 


1.94 


38.7 


1.07 


19.95 


551.55 


1.025 


10 9 GeV 


2.28 


45.3 


1.60 


19.87 


701.75 


1.062 


2 x 10 16 GeV 


1.33 


26.5 


1.00 


19.92 


751.88 


1.072 



interesting phenomenon, and indicates that Koide's relation is a universal result. 

Furthermore, we can find from Table II that k u is not 1 as ki of charged leptons. This 
means that Koide's relation should be improved before being applied to the case of quarks. 
In we calculated the range of k u , and got 1.1 < k u < 1.4 (with the mean value o 1.25). 
However, in we did not consider the renormalization effect of quark masses. With this 
effect taken into account, we find that k u changes a little, from 1.25 to 1.33. 

From Table III, we can find that kd also keeps invariant with the change of energy scales, 
just as k u , and its approximate value is 1.05. 

Similarly, in Table IV, k\ with the increase of /i at high energy scales are listed, both 
in the standard model and in the minimal SUSY model, as we know that the values in ?J 
should only be taken at low energy scale. 

At the same time, we can see in Table IV that ki at high energy scales is still quite close 
to 1 as at the low energy scale, which means that Koide's relation is suitable for charged 
leptons at all the energy scales, just as quarks. Also, we can find in Table IV that k\ in the 
minimal SUSY model is a little larger than that in the standard model. 
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TABLE IV: xi, y\ and k% at different energy scales. The upper three rows are the cases in the 
standard model, and the lower three rows are the cases in the minimal SUSY model. 



A 1 


m e (MeV) 


m M (MeV) 


m T (GeV) 




Vl 


h 


91 GeV 


0.48684727 


102.75138 


1.7467 


211.05465 


3587.78 


1.001881 


10 9 GeV 


0.51541746 


108.78126 


1.8492 


211.05467 


3587.77 


1.001881 


2 x 10 16 GeV 


0.49348567 


104.15246 


1.7706 


211.05468 


3587.95 


1.001888 


91 GeV 


0.48684727 


102.75138 


1.7467 


211.05465 


3587.78 


1.001881 


10 9 GeV 


0.40850306 


86.21727 


1.4695 


211.05661 


3597.28 


1.002277 


2 x 10 16 GeV 


0.32502032 


68.59813 


1.1714 


211.05797 


3604.08 


1.002560 



In summary, we can conclude that Koide's relation is a result independent (or insensitive) 
of energy scales of interaction, from low to extremely high energies. 

III. KOIDE'S RELATION FOR NEUTRINOS 

After testing Koide's relation for quarks and charged leptons, and finding that k u , 
and ki are independent of energy scales, a natural question is what about this relation 
for neutrinos. Here we again introduce the parameter k v as in Ref. Ill ], and discuss this 
problem. Moreover, we try to find the relations between the four parameters k u , kd, k v and 
ki, and finally we can get the neutrino masses with some conjectures. 

It is quite difficult to verify whether Koide's relation is suitable for neutrinos because of 
the long-existed inaccuracy of the experimental data of neutrinos. Due to the untiring efforts 
by the numerous neutrino experiments, the oscillations and mixings of neutrinos have been 
strongly established now. The solar neutrino deficit is caused by the oscillation from v e to a 
mixture o f v v and v T with a mixing angle approximately of 9 so \ « 34° in the KamLAND 
and SNO [lq experiments. Also, the atmospheric neutrino anomaly is due to the to v T 
oscillation with almost the largest mixing angle of a t m ~ 45° in the K2K 17] and Super- 
experiments. However, the non-observation of the disappearance of v e in 
1£ ] experiment showed that the mixing angle 6^ is smaller than 5° at the best 



Kamiokande 
the CHOOZ 
fit point 



21 



S 



These experiments also measured the mass-squared differences of the neutrino mass eigen- 
states. Accor ding to the global analysis of the experimental results, we have (the allowed 



ranges at 3a) 



21] 



1.4 x 1(T 3 eV 2 < Am 2 tm = \m\ - m\\ < 3.7 x 1(T 3 eV 2 , (5) 



and 

5.4 x 1(T 5 eV 2 < Am 2 ol = \m\ - m\\ < 9.5 x 1(T 5 eV 2 , (6) 

where m 1; m 2 , 777.3 are the masses of the three mass eigenstates of neutrinos, and the best 
fit points are |m§ - m\\ = 2.6 x 1(T 3 eV 2 , and \m\ - m\\ = 6.9 x 10" 5 eV 2 |2l|. 

Due to Mikheyev-Smirnov-Wolfenstein |22j matter effect of solar neutrinos, we already 



know that 77i 2 > m\. Hence we have 

m\ = m 2 2 - Am s 2 ol , (7) 

and 

m\ = m\± Am 2 tm . (8) 

So there are two mass schemes, (1) the normal mass scheme 777,3 > ^2 > mi, and (2) the 
inverted mass scheme 7772 > mi > 7773. We will discuss both of them in the following. 

Now we apply Koide's relation to neutrinos. First we take the normal mass scheme for 
example. If Koide's relation holds well for neutrinos, we have 

2 

mi + m 2 + m 3 = - ( v /m7+ -Jm^+ y/mj) . (9) 
o 

Substituting Eqs. © and © into Eq. ©, we get 



m 2 2 - Am 2 sol + 777 2 + J ml + A7?7 2 tm = - {/ m\ - Am s 2 ol + ^/m~ 2 +t/m% + Aml tm 



Solving this equation, we find that there is no real root for m 2 with the restrictions in 
Eqs. (0) and ©. This means that no matter what value 777 2 is, Koide's relation does not 
hold for neutrinos. So is the inverted mass scheme. 

Thus we must introduce a parameter k u to character the deviation of neutrinos from 



Koide's relation, 



mi + m 2 + m 3 

kv = 2/ / T / T / A 2 - ( 10 ) 
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FIG. 1: The range of k u of the normal mass scheme 7713 > m<i > m\. 

Since k u ^ 1, we must determine its range, and this can help us to find the relations 
between k u , kd, k v and ki, and fix the neutrino masses. We now check the situations for the 
two mass schemes, respectively. 

1. For the normal mass scheme, m 3 > m 2 > m 1; we have 



k = ~ A^soi + m 2 + yjml + Am^ tm 

I (^-Am s 2 ol + v^+ \/«2 + A^atm) 

We can see that k v is the function of m 2 if Am s 2 ol and Am 2 tm are fixed. Due to the inaccuracy 
of the experimental data, we take Am 2 ol and Am 2 tB as their best fit values here. The range 
of k v is shown in Fig. 1. 

We can see from Fig. 1 that 0.50 < k v < 0.85, and k v decreases with the increase of m?,- 
So k v < 1 for neutrinos. 

2. For the inverted mass scheme, m 2 > mi > 777,3, we have 



k = Vmj- Am 2 sol + m 2 + ^m\- Am 2 tm ^ 

I (\/™2- Aml, + \/^+ V /m 2 - Am atm) 

The range of A;^ is shown in Fig. 2. 

We can see from Fig. 2 that 0.50 < k v < 0.65. Again, k„ < 1 for neutrinos. 

Altogether, 0.50 < k u < 0.85 for both of these two mass schemes. And k v of the normal 
mass scheme is larger than that of the inverted mass scheme. 
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FIG. 2: The range of k u of the inverted mass scheme mi > rrt\ > 7713. 
Conclusively, the values of k u , kd, k v and k\ can be summarized as follows, 

and MMf«U» 

\dj \s J \bj k d ^l 

We believe that the problem of the origin of the masses of leptons must be solved together 
with that of quarks. Since ki — 1 and kd « 1, we may conjecture that k\ + kd ~ 2. At the 
same time, since 0.50 < k u < 0.85 and fc u ~ 1.33, we may analogize the conjecture of fc/ and 
fed, and propose the hypothesis that 

K + k u ^ 2. (14) 

This is from the speculation that there must be some relation between k u , kd, k v and k\. 
Of course, this Ansatz is not the unique one of the relations between k u , kd, k v and k\. For 
example, we may also assume that 

kik d w k v k u w 1, (15) 

or 

k\ kd k v k u 

Eq. (fT^j) is based on the assumption that 0i + 9 d ~ U + 9 U ?z | in Foot's geometrical 
interpretation. 

All these hypotheses can show the balance between fc„ and fc u intuitively and trans- 
parently. This situation seems to be similar to the quark-lepton complementarity between 
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mixing angles of quarks and leptons 23|, and we may call it a quark- lepton complementarity 
of masses. 

From Table II, we can see that the value of k u is approximately 1.33. Thus from the 
hypothesis k v + k u « 2, we get that k v w 0.67. This is consistent with the normal mass 
scheme and in conflict with the inverted mass scheme. This indicates that the three masses 
of neutrino mass eigenstates are heavier in order, which is the same as quarks and charged 
leptons. 

Now we can estimate the absolute masses of neutrinos. Substituting k u = 0.67, Am^ tm = 
2.6 x 10~ 3 eV 2 , and Am 2 ol = 6.9 x 10~ 5 eV 2 into Eq. (fTTjl . we obtain the value of m 2 , 



og7 _ y/m% - 6.9 x 10" 5 eV 2 + m 2 + ^Jm\ + 2.6 x 10" 3 eV 2 

§ [\Jm\ - 6.9 x 10" 5 eV 2 + jm~ 2 + tym\ + 2.6 x 10" 3 eV 2 ) 2 



and we get m 2 = 0.0089 eV. 
Straightforwardly, we get 



77H = J m 2 2 - Am s 2 ol = 0.0032 eV, (18) 

and 



m 3 = ^J m \ + Aml tm = 0.052 eV. (19) 

Now we discuss the uncertainty of mi, m 2 and 1713. In Fig. 1, we can see that the slope 
of the curve in very large where 0.65 < k v < 0.8, so the value of m 2 is not very sensitive 
to the error of k v . m 2 will approximately be 8.5 ~ 8.9 x 10~ 3 eV even if the value of k v 
charges from 0.65 to 0.8, so the value of m 2 is precise to a fairly good degree of accuracy. 
Similarly, the value of m 3 will be about 0.052 eV to a good degree of accuracy too, because 
m 3 = y/m\ + Am^ tm , and Am^ tm ^ m|. The only point desired to be mentioned here is 
the range of mi. Because if m\ is rather closed to Am 2 ol , and due to the big uncertainty 
of Am 2 ol , the value of mi may change largely with k v . We can see this in the other two 
hypotheses in Eqs. (fT3j) and (JTfij) . 

1. In Eq. (fT5|). where k\kd ~ k v k u w 1, we have k v w 0.75, and 

mi = 0.0012 eV, 
m 2 = 0.0084 eV, 

m 3 = 0.050 eV. (20) 



12 



2. In Eq. (ITBl . where ^- + Tr~r- + -r-~2, we have k v m 0.80, and 

mi = 1.0 x 10" 5 eV, 
m 2 = 0.0084 eV, 

m 3 = 0.050 eV. (21) 

From Eqs. (fTSj) . (f2~Uj) and (}2~T]) . we can see that the value m x is only a rough estimate of 
the first step till now, and its effective number and order of magnitude may change with the 
more and more precise experimental data in the future. However, the values of m 2 and m% 
are consistent in these three hypotheses. 

Finally, we should point out that Koide |24] also gave an interpretation of his relation as 
a mixing between octet and singlet components in a nonet scheme of the flavor U(3), and 
got the neutrino masses as mi = 0.0026 eV, m 2 = 0.0075 eV and m^ = 0.050 eV [25j. We 
can see that his results are strongly consistent with ours, especially with the results from 
Eq. (I14j) . And recently, he got the three masses as mi = 0.0039 eV, m 2 = 0.0088 eV and 
m 3 = 0.053 eV in a seesaw mass matrix model of quarks and leptons with flavor-triplet 



Higgs scalars 



which is even closer to our results from Eq. (fTUh 



IV. SUMMARY 

Finally, we give some discussions. 

1. Since Koide's relation is such a wonderful result for charged leptons at low energy 
scale, to explore its behavior at high energy is straightforward. We carefully test whether it 
is energy scale independent in Section II, and find that it is really independent or insensitive 
of energy scale in a huge energy range and is almost the same in both the standard model 
and the minimal SUSY model, which proves that Koide's relation is a universal result in 
particle physics. 



2. Seesaw mechanism 



271 ] may give the origin of neutrino masses, in which the right- 



handed very- heavy neutrinos is included into the Lagrangian of the standard model, i.e., 




(22) 
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where the scale of mo is characterized to be the energy scale v of the electroweak spontaneous 
breaking, and Mr is the right-handed very-heavy neutrino mass matrix. 
From Eq. (J22j) . we can get 

m v = -m D M R mJ), 

where the neutrino mass eigenvalues are mi = m„/Mi, m 2 = m^jM-i and = m^/M^. 
Thus, the smallness of neutrino masses is due to the large values of M%, M 2 and M 3 . (For 
example, m u ~ 0.1 eV if Mr is taken as 10 14 GeV.) 

However, seesaw mechanism can only present an illustrative interpretation of the origin 
of neutrino masses, without accurate prediction of the masses of neutrino mass eigenstates. 
To obtain those, we must extend our theory and make some speculation, i.e., we examine 
whether Koide's relation of charged leptons also holds well for neutrinos, and we find that 
not all quarks and leptons obey Koide's relation precisely. So we introduce the parameters 
k u , kd and k v to describe the deviations of quarks and neutrinos from the exact Koide's 
relation. With this improved relation and the conjecture of a quark-lepton complementarity 
of masses such as ki + kd ~ k v + k u ~ 2, kikd ~ k v k u ^lor- 1 --|--^-~r--|-r-~2, we 
can determine the absolute masses of the neutrino mass eigenstates. Due to the inaccuracy 
of the experimental data of neutrinos nowadays, these results (especially the value of mi) 
should be only taken as primary estimates. It is also possible that seesaw mechanism is 
responsible for the deviation of the neutrino masses from the exact Koide's relation. 

3. There remain some open questions to be answered. Such as 

(1) . Which hypothesis between k u , kd, k v and ki is really the relation between them? 
For example, if^- + r-~r- + r-~2, we have Q\ + 6 d ~ 6 U + U « £ in Foot's geomet- 

KI Kd Ku K u 2 

rical interpretation, is there really some deeper reason behind it, just like the quark-lepton 
complementarity of their mixing angles 

(2) . We can see from Table II that k u is approximately 1.33, so may it be | exactly? If 
so, we have k v = | in Eq. (JHJ), ^ = f m Eq. (JUJ), and k u = | in Eq. (fTBj) . All these k u 
and k v are of special values, and is there some unknown principle leading this? Moreover, 
is kd — 1 exactly as ki, or deviates from 1 slightly? If so, why? 

In conclusion, we believe that there must be some deeper principle behind the elegant 
Koide's relation for charged leptons, and thus it is meaningful to check whether this relation 
is also applicable to quarks and neutrinos, at least at some degree. We show in this paper 
that the Koide's relation with its deviation characterized by the parameters k u and kd is 
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also applicable to quarks without sensitivity to energy scale. By using the improved Koide's 
relation with an Ansatz of quark-lepton complementarity of masses, we can also determine 
neutrino masses. If the prediction of neutrino masses is tested to be consistent with the 
precise experiments in the future, it would be a big success of Koide's relation, which may 
shed light on our way to a grand unification theory of quarks and leptons. 
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